Draft version April 6, 2012 

Preprint typeset using I^T^^X style emulatcapj v. 5/2/11 



SN2010AY IS A LUMINOUS AND BROAD-LINED TYPE IC SUPERNOVA WITHIN A LOW-METALLICITY 

HOST GALAXY 

N. E. Sanders\ a. M. Soderberg\ S. Valenti^, R. J. Foley\ R. Chornock\ L. Chomiuk^-^, E. Berger\ S. 
Smartt^ K. Hurley^, S. D. Barthelmy^ E. M. Levesque'', G. Narayan"^, M.T. Botticella^ M. S. Briggs*, V. 
Connaughton'*, Y. Terada'^ N. Gehrels'', S. Golenetskii"', E. Mazets^", T. Cline", A. von Kienlin^^, W. 
Boynton^-', K. C. Chambers", T. GRAv^^ J. N. Heasley", K. W. Hodapp", R. Jedicke", N. Kaiser", R. P. 
Kirshner\ R.-P. Kudritzki", G. A. Luppino", R. H. Lupton^'', E. A. Magnier", D. G. Monet", J. S. Morgan", 
P. M. Onaka", p. a. PRICE^^ C. W. Stubbs^, J. L. Tonry", R. J. Wainscoat", M. F. Waterson^*. 

Draft version April 6, 2012 

ABSTRACT 

We report on our serendipitous pre-discovery detection and detailed follow-up of the broad-lined 
Type Ic supernova (SN) 2010ay at z = 0.067 imaged by the Pan-STARRSl Stt survey just 4 
days after explosion. The SN had a peak luminosity, Mr « —20.2 mag, significantly more luminous 
than known GRB-SNe and one of the most luminous SNe Ib/c ever discovered. The absorption 
velocity of SN2010ay is vsi w 19 x 10'^ km s~^ at ^ 40 days after explosion, 2 — 5 times higher 
than other broad-lined SNe and similar to the GRB-SN2010bh at comparable epochs. Moreover, 
the velocity declines ^ 2 times slower than other SNe Ic-BL and GRB-SNe. Assuming that the 
optical emission is powered by radioactive decay, the peak magnitude implies the synthesis of an 
unusually large mass of ^^Ni, A/ni = 0.9 Mq. Modeling of the light-curve points to a total ejecta 
mass, Afoj ~ A.7Mq, and total kinetic energy, Ek ~ 11 x 10'"'^ ergs. The ratio of A/ni to Afcj is ~ 2 
times as large for SN 2010ay as typical GRB-SNe and may suggest an additional energy reservoir. The 
metallicity (Iog(0/H)ppo4 -I- 12 = 8.19) of the explosion site within the host galaxy places SN2010ay 
in the low-metallicity regime populated by GRB-SNe, and ^ 0.5(0.2) dex lower than that typically 
measured for the host environments of normal (broad-lined) Ic supernovac. We constrain any gamma- 
ray emission with <6 x 10"** erg (25-150 kcV) and our deep radio follow-up observations with the 
Expanded Very Large Array rule out relativistic ejecta with energy, E > 10*® erg. We therefore rule 
out the association of a relativistic outflow like those which accompanied SN 1998bw and traditional 
long-duration GRBs, but place less-stringent constraints on a weak afterglow like that seen from 
XRF 060218. These observations challenge the importance of progenitor metallicity for the production 
of a GRB, and suggest that other parameters also play a key role. 

Subject headings: Surveys:Pan-STARRSl — gamma-rays: bursts — supernovac: individual (2010ay) 
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1. INTRODUCTION 

Recent observations have shown that long-duration 
gamma-ray bursts are accompanied by Type Ic super- 
novae (SNe) with broad absorption features (hereafter, 
"broad- Uned," BL), in dicative of high photosp heric ex- 
pansion velocities (see lWooslev fc Bloom|[2006l for a re- 
view). This GRB-SN connection is p opularly explained 
by t he favored "collapsar model" (jMacFadven et aH 
|2001| ) in which the gravitational collapse of a massive 
(M > 20 Mq) progenitor star gives birth to a central 
engine - a rapidly rotating and accreting compact object 
- that powers a relativistic outflow. At the same time, 
not all SNe Ic-BL show evidence for a central engine. 
Radio observations constrain the fraction of SNe Ic-BL 
harboring relativistic outflow s to be less than a third 
(jSoderberg et al.l[20?)6al IMH . 

The physical parameter(s) that distinguish the pro- 
genitors of GRB-associated SNe from other SNe Ic-BL 
remains debated, while theoretical considerations indi- 
cate that progenitor metallicity may play a key role 
(jWooslev fc I Ieger.,2006i) . In the collapsar scenario, mas- 
sive progenitor stars with metallicity above a thresh- 
old, Z > 0.3 Zq, lose angular momentum to metal 
line-driven winds, preventing the formation of a rapidly 
rotating compact remnant, and in turn, a relativistic 
outflow. At the same time, the hydrogen-free spec- 
tra of SNe Ic-BL indicate that their stellar envelopes 
have been stripped prior to explosion, requiring higher 
metalli cities (e.g., Z » Zc^ ) if due to radiation driven- 
winds (jWooslev et al.|[T995D . Alternatively, short-period 
(~ 0.1 days) binary interaction may be invoked to spin up 
stars via tidal forces as well as ca use mass loss via Roche 
lobe overflow (jPodsiadlowski et al . 2004: Frver fc Hegerl 
I2005f ). However, even in the binary scenario, GRB for- 
mation is predicted to occur at higher rates in lower- 
metallicity environments, where the radius and mas s loss 
rates of stars should be smaller ([Izzard et al.l 120041 ). 

Observationally, most GR B-SNe are discovered within 
dwarf star- forming galaxies (jFruchter et al.l [20061 ) char- 
acterized by sub-sola r metallicities, Z < 0.5 Zq 
(|Levesaue et al.|[2"010al ) . This has been interpreted as ob- 
servational support for the me tallicity-dependent collap- 
sar model (jStanek et al.ll2006D . Meanwhile, SNe Ic-BL 
without associated GRBs hav e historically been f ound in 
more enriched environments (jModjaz et al.l[2008h . moti- 
vating the suggestion of an observationally determined 
"cut-off metallicity" a bove which GRB-SNe do not form 
(jKocevski et al.ir2009[ ). However, this difference may be 
partly attributable to the different survey strategies: SNe 
have been found in large numbers by galaxy-targeted sur- 
veys biased towards more luminous (and therefore higher 
metallicity) environments, while GRB host galaxies are 
found in an untargeted manner through their gamma-ray 
emission. 

Against this backdrop of progress, recent observations 
have begun to call into question some aspects of this 
scenario. First, several long-duration GRBs have now 
been identified in solar or s uper-solar metallicity e nviron- 
ments (e.g., GRB 020819: iLevesaue eralll2010b[) . Sim- 
ilarly, the luminous radio emission seen from SN Ic- 
BL 2009bb pointed unequivocally to the production 
of co pious relativistic ejecta resembling a GRB after- 
glow (jSoderberg et al.l 120101 ) while the explosion envi- 



ronment was characterized by a sup er-solar metallicity, 
Z - 1 - 2 Zq (jLevesque et alboiOd) . Together with the 
growing lack of evidence for massive progenitor stars for 
SNe Ic in pre -explosion Hubble Space Telescope images 
(|Smarttll200"9f ). a lower mass (M ~ 10 - 20 Mq) binary 
progenitor system model (with a gentler me tallicity de- 
pend ence) is gaining increasing popularity (|Yoon et al.l 
|2010[ ). Multi- wavelength studies of SNe Ic-BL in metal- 
poor environments may shed further light on the role of 
metallicity in the nature of the progenitor and the explo- 
sion properties, including the production of relativistic 
ejecta. 

Fortunately, with the recent advent of wide-field opti- 
cal transient sur veys (e.g., Catalin a Real-time Transient 
Survey fCRTS): [Drake et al.|[200l Panchromatic Survey 
Telescope and Rapi d Response System ; (Pan-STARRSl, 
abbreviated P SD iKaiser et aD 120021 . Palomar Tran- 
sient Factory; iLaw et al .1 I2009D SNe Ic-BL are now 
being d iscovered in meta l -poor environments , Z ~ 
0.5 Zq (jArcavi et al.l 120101: iModiaz et aI[|20lH) thanks 
to an unbiased search technique. In this paper, we 
present pre-discovery Pan-STARRSl imaging and multi- 
wavelength follow-up observations for the broad-lined 
Type Ic SN2010ay discovered by CRTS ([Drake et all 
I2010f) . In 31 we report our optical (Pan-STARRSl, 
Gemini, William Herschel Telescope) and radio (Ex- 
panded Very Large Array; EVLA) observations. In ij3l 
we model the optical light-curve and analyze the spec- 
tra to derive the explosion properties of SN2010ay. In 
Sin we use our observations of SN 2010ay with the EVLA 
to place strict limits on the presence of relativistic out- 
fiow. In ij5l we draw from gamma-ray satellite coverage 
to rule out a detected gamma-ray burst in association 
with SN 2010ay. In fjGl we derive the explosion site metal- 
licity and find it to be significantly sub-solar and typical 
of most GRB-SNe host environments. In ^J7l we discuss 
the implications of these findings in the context of the 
explosion and progenitor properties and conclude in SjSj 

2. OBSERVATIONS 
2.1. Discovery by CRTS 

SN2010ay was discov ered by the Cata lina Real- 
time Transient Su rvev ([Drake et al I I2009D on 2010 
March 17.38 UT ([Drake et al.l I2010D and designated 
CSS100317:123527-H270403, with an unfiltered magni- 
tude of rn w 17.5 mag and located < 1" of the cen- 
ter of a compact galaxy, SDSS J123527. 19+270402.7 at 
z = 0.067 (Table tt]). We adopt a distance, Dl ^ 297.9 
Mpc, to the host galax}{3, and note that the Galac- 
tic extinction toward this galaxy is E{B — V) ~ 0.017 
(jSchlegel et al.l[T998h . Pre-d iscovery unfiltered images 
from CRTS revealed an earlier detection of the SN on 
Mar 5.45 UT at m « 18.2 mag a nd a non-detection from 
Feb 17.45 UT at m > 18.3 mag ([Drake et al.l l2010'). 

A spectrum obtained on Mar 22 UT revealed the 
SN to be of Type Ic with broad features, similar to 
the GRB-associa ted SN 1998bw spectru m obtained near 
maximum fight (iFflippenko et all I2010D . This classifi- 
cation was confirmed bv lPrietol ( 20101 ) who additionally 
reported photometry for the SN (see Table [2|). After 
numerically subtracting the host galaxy emission, they 

I'J We assume Ho = 71 km s^^ Mpc^i, Qa = 0.73, Qm = 0.27 
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TABLE 1 

SN 2010AY HOST GALAXY SDSS 

J123527. 19+270402. 7 



Parameter 


Value 


RA 


12h35'"27M9 (J2000) 


DEC 


+27°04'02.7" (J2000) 


Redshift (z) 


0.0671 


Petrosian radius 1.355" 


Photometry'' 


u' 


19.56 ± 0.03 mag 


g' 


19.02 ± 0.01 mag 


r' 


19.02 ± 0.01 mag 


i' 


18.69 ± 0.01 mag 


z' 


18.87 ± 0.04 mag 


U 


19.50 ±0.06 mag 


B 


19.02 ± 0.05 mag 


V 


19.13 ± 0.05 mag 


R 


18.94 ± 0.05 mag 


I 


18.90 ± 0.06 mag 


Extinction 


E(B-V)mw'' 


0.017 mag 


E(B-V)hoBt" 


0.2 mag 



Note. — SDSS host galaxy properties and 
ugriz photometry. 

^ Model magnitu des from SDSS DR8 
JAihara et al.| 120111) . Host galaxy photom- 
etry has not been corrected for extinction. 
U BY RI photometry has been converted from 
the SDSS ugriz measurements using the 
transformation of lBlanton RoweisI < |2007D . 
^ The Milky Way extinction as determined by 
ISchlegel eraTI {1991), assuming Rv — 3.1. 
^ The host galaxy extinction determined from 
the SDSS spectrum centered on the galaxy 
nucleus, via the Balmer decrement as de- 
scribed in !|2.3.1I 

estimate an unusually luminius absolute magnitude of 
V « -19.4 mag. 

2.2. Pre-Discovery Detection With Pan-STARRSl Stt 

The field of SN2010ay was screndipitously observed 
with the PSl Stt survey in the weeks preceding its 
discovery. Pan-STARRSl is a wide-field imaging sys- 
tem at Haleakala, Hawaii dedicated to survey obser- 
yations (iKaiser et al.l |2002[ ). The PSl optical design 
(|Hodapp et all I2004D uses a 1.8 meter diameter //4.4 
primary mirror, and a 0.9 m secondary. The telescope il- 
luminat es a diameter of 3.3 deg rees. The Pan-STARRSl 
imager (|Tonrv fc Onakal [2001 comprises a total of 60 
4800 X 4800 pixel detectors, with 10 fim pixels that sub- 
tend 0.258 arcsec. The PSl observations are obtained 
through a set of five broadband filters designated as gpi, 
''PI, *pi, -Zpi, and ypi. These filters are similar to those 
used in previous surveys, such as SDSS (jFukugita et al.l 
|1996[) . However, The ^pi filter extends 20 nm redward 
of 5', the zpi filter is cut o ff at 930 nm, and S PSS has 
no corresponding ypi filter (jStubbs et al.ll2010l ). 

The field of SN2010ay was observed on 2010 Febru- 
ary 21st (rpi-band) and February 25th (ipi-band. Fig- 
ure [T]). On each night, four exposures were collected fol- 
lowing the strategy of the PSl-3pi survey (Chambers et 
al., in preparation). Following the CRTS discovery and 
announcement of SN2010ay, we geometrically registered 
SDSS pre-explosion images to the PSl images and per- 



formed digital image subtraction using the ISIS package 
(|Alardll2000D . No residual flux was found in the differ- 
ence rpi-band image from 2010 February 21 with an up- 
per limit of r > 22.0 mag. However, we detect residual 
flux at the position of SN2010ay in the ipi-band resid- 
ual image from 2010 February 25 with a magnitude of 
i = 21.1 ±0.3 mag. 

The field was again observed in the ipiband on 2011 
March 25 and the rpi-band on 2011 March 29, but no 
residual flux was detected in the subtractions at the SN 
position to limits of i' > 22.2 and / > 21.9. 

In Table [2] and Figure [21 we compile photometry from 
the PSl detections, our optical observations, and the cir- 
culars to construct a light-curve for SN2010ay. 

2.3. Optical observations 

We obtained an optical spectrum ('^ 3000—11000 A) of 
SN 2010ay on April 1 UT, from the ISIS blue arm instru- 
ment of the 4.2 m Wflham Herschel Telescope (WHT) 
at the Roque de Los Muchachos Observatory. The spec- 
trum was taken at the parallactic angle and the expo- 
sure time was 1800 sec. We obtained a second 1800 sec 
optical spectrum (~ 3600 — 9600 A) using the Gemini 
Muhi-Object Spectrograph (CMOS) on the 8.1 m Gem- 
ini North telescope on 2010 April 11.4 UT. We employed 
standard two-dimensional long-slit image reduction and 
spectral extraction routines in IRAlH- We do not apply 
a correction for atmospheric differential refraction, be- 
cause the displacement should be < 0.5" at the airmass 
of the observations, « 1.0. 

In both our Gemini and WHT spectra, broad absorp- 
tion features associated with the SN are clearly detected 
in addition to narrow emission lines typical of star form- 
ing galaxies. We distinguish the host galaxy emission 
from the continuum dominated by the highly broadened 
SN emission by subtracting a high-order spline flt to 
the continuum. Both SN and host galaxy spectral com- 
ponents are shown in Figure [3l As illustrated in the 
Figure, the broad, highly-blended spectral features of 
SN 2010ay resemble those of the Type Ic-BL SN 2010bh 
associated with GRB1 00316D) at a similar similar epoch 
Chornock et al.l [2010D . In particular, the broadening 
and blueshift of the feature near 6355 A are similar for 
SN2010ay and SN2010bh, and are broader and more 
blueshifted than in SN 1998bw at a comparable epoch. 
We discuss the comparison between these two SN fur- 
ther in gSJ and 

Additionally, we obtained 60 s gri-hand images of 
SN2010ay on 2010 April 11.4 UT using CMOS. The 
data were reduced using the gemini package in IRAF, 
and photometry was performed using the standard 
CMOS zero-point^. We measure that [g,r,i] = 
[18.90, 18.32, 18.04] ±0.1 mag. 

Imaging photometry is not available at the epoch 
of our WHT spectrum. However, the spectrum was 
flux-calibrated against observations of the standard star 
Feige34, which was observed the same night and at ap- 
proximately the same low airmass as the supernova. For 

'^^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

http://www. gemini. edu/?q=node/10445 
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iModiaz et all (|20Tol) was also similar: E(B-V)=: 0.3 mag. 
Furthermore, we note that the co lo r (B — V = 0.78 ± 
0.07 mag) as reported by IPrietol (|2010l ) at 4 d after 
i?— band peak is significantly redde r than SN Ib/c color 
curve templates iDrout et al.l ()2011|) , further supporting 
a non-negligible host galaxy extinction. 

2.4. Radio Observations 



Fig. 1. — Images illustrating the Pa n-STARRSl prc-disco vcry detection of SN2010ay and the surrounding field, a) Pre-explosion gri- 
composite image from the SDSS DR7 (lAbazaiian et al.ll2009l ). observed 2004 December 21, b) SDSS i-band image geometrically registered 
to the PSl image frames and zoomed in to the host galaxy of SN 2010ay, c) ipj-band image from the Stt survey of PSl, observed 2010 
February 25, d) the difference of the SDSS i' and PSl ipi images. The transient omission can be seen in the NE corner of the last panel. 
Nearby stars are included in the zoomed-in frames to illustrate the efficiency of subtraction. 

the observations of both the standard star and super- 
nova, the slit was placed at the parallactic angle. We 
perform synthetic photometry on the spectrum to extract 
the flux at the central frequency of the R-band (6527 A) 
and find R = f 8.2 ± 0.2 mag. We then subtract the host 
galaxy flux numerically. 

2.3.1. Host galaxy features 

We measure fluxes of the narrow emission lines from 
the host galaxy in our Gemini spectrum, as reported in 
Table [31 We fit a Gaussian profile to each narrow emis- 
sion line; for nearby lines such as [N II] and Ha, we fit 
multiple Gaussians simultaneously. We model the local 
continuum with a linear fit to 20 A regions on either side 
of each line. We estimate uncertainties in quantities de- 
rived from the line fluxes by Monte Carlo propagation of 
the uncertainties in the flux measurement. 

The host galaxy is signiflcantly reddened as evidenced 
by the flux ratio of H„ to Bp, w 3.45 ± 0.02. We infer 
E(B — V) = 0.2 mag {Ay = 0.6 mag), as measured from 
the Balmer decrement in our Gem ini spectrum, assuming 
Ry = 3.1, Case B recombination ( Qsterbrock fc Ferlandl 
I2006f) . and the reddening law of lCardelh et al.l (|1989D 
This is similar to the value derived from the SDSS 
DR8 nuclear fiber spectroscopy line fiuxcs for the host 
galaxy (E(B — V) = 0.2 mag). The value reported in 



We observed SN2010ay with the EVLA (jPerlev et al.l 
12001 on three epochs, 2010 March 26, 2010 April 29, 
and 2011 May 7. All EVLA observations were obtained 
with a bandwidth of 256 MHz centered at 4.9 GHz. We 
used calibrator J1221+2813 to monitor the phase and 
3C286 for fiux calibration. Data were reduced using the 
standard packages of the Astronomical Image Processing 
System (AlPS). We do not detect a radio counterpart to 
SN2010ay in these observations and place upper limits 
of Fi, < 46,42,30 /xJy (3(t) for each epoch respectively 
corresponding to upper limits on the spectral luminosity 
spanning i,y < (3.6 — 5.5) x 10^^ erg cm~^ s~^. 

As shown in Figure SI these limits are compa- 
rable to the p eak luminosities observed for ordi- 
nary SNe Ib/c (iBerger et all I2003at ISoderberd [20071: 



. erge . . _ . . . . 

[Chevalier fc FranssonI 120061: iSoderberg et al.l 120101 and 
references within) and a factor of 10^ to 10^ less luminous 
than the radio afterglows associated with GRBs 020903, 
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TABLE 2 

SN 2010AY LIGHT-CURVE PHOTOMETRY 



UT date 


MJD 


J. a 


Filter 






Source 


733820.4 


55244.4 


-29 




< 18.3 


< -19.2 ±0.2 


d 


733824.6 


55248.6 


-25 




< 22.0 ±0.1 


< -16.0 ±0.1 


e 


733828.2 


55252.2 


-21 




21.1 ±0.3 


-16.8 ±0.3 


e 


733830.4 


55260.4 


-13 




18.2 


— 19.3 ± 0.2 


d 


733848.4 


55272.4 


-1 




17.5 


-20.4 ±0.2 


d 


733853.2 


55277.2 


4 


B 


18.39 ±0.05 




f 


733853.2 


55277.2 


4 


V 


17.61 ±0.05 




f 


733853.2 


55277.2 


4 


R 


17.44 ± 0.05 


-20.22 ±0.07 


f 


733863.0 


55287.0 


14 


R 


18.2 ±0.2 


-19.2 ±0.2 


g 


733873.0 


55297.0 


24 


g 


18.9 ±0.1 




h 


733873.0 


55297.0 


24 


r 


18.3 ±0.1 


-19.0 ±0.1 


h 


733873.0 


55297.0 


24 


i 


18.0 ±0.1 




h 


734221.0 


55645.0 


372 




< 22.2 ±0.2 


< -16.22 ±0.10 


c 


734225.0 


55649.0 


376 


rpi 


< 21.9 ±0.1 


< -15.7±0.1 


c 



^ Time since peak in days, relative to the fitted value: 2010 March 18.5 it 0.2. 

^ The measured apparent magnitude of the source, in the filter noted and without extinction 
correction. For the Pan-STARRSl photometry, a template image was subtracted; for the other 
points, the host galaxy fiux has not been subtracted. 

^ The absolute R magnitude of the SN. Filter conversion, host flux subtraction, and extinction 
correction have been performed as described in ^3.11 
^ IDrake eraTI ( I2010D . unfiltered (synthetic V-band). 
^ Photometry from the Pan-STARRSl Stt survey. 
^ iPrietol II2010D . 

^ Synthetic photometry obtained from our WHT spectrum as described in ^2.31 
^ Photometry from our Gemini/GMOS observations described in ^2.2| 
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Fig. 2. — The optical R-band light curve of SN2010ay, as com- 
piled in Table|2]from CBET 2224 (red squares), the PSl Stt s urvey 
(blue), our Gemini images (black), and synthetic photometry ( i|2.3l l 
based on our WHT spectrum (gold). Triangles denote upper limits. 
The thick dashed line represents the lumi nosit y of an expanding 
fireball fit to our early-time photometry ( i|3.1ll. The thin da shed 
line is the SN Ib/c light-curve template of lDrout et al.l 1 1201 II ) and 
the gray field represents the standard deviation among its con- 
stituent photometry. The template is stretched by (1 + 2:) = 1.067 
with the best fit parameters tRpeak = Marchl8±2 (MJD 55273±2), 



Rpcak 



-20.2 ± 0.2 mag. The solid line is the light-curve of 



the GRB-SN 1998bw IIGalama et aLlirggg ). redshifted to match 
SN 2010ay. The vertical dotted lines mark the epochs of our Gem- 
ini and WHT spectroscopy. 



030329, and 031203 at ea rly epochs (iBerger et al.l l2003bl: 
ISoderberg et al.ll2004allbl : iFrail et alll2005D . In compar- 
ison with the radio l uminosities observed for the rcla- 
tivistic SNe 1998bw (p<:ulkarni et al.l I1998D and 2009bb 
(jSoderberg et al.ll2010[) . SN2010avis a factor of > 10 less 



TABLE 3 

Emission line fluxes measured for the 
host galaxy of sn2010ay 



Emission Line 


Flux 




^]^Q— I6gj.g s~-'-cm~^) 



[O II]A3726,3729 

H5 

H7 

[O III]A4363 
H/3 

[O III]A4959 
[O III]A5007 
[N II]A6548 
Hq 

[N H]A6584 
[S H]A6717 
[S H]A6731 



52 ± 1 
4.5 ±0.2 
9.9 ±0.1 
0.7 ±0.1 
24.5 ±0.1 

30.4 ±0.1 
90.1 ±0.2 

2.2 ±0.1 

84.5 ±0.1 
6.33 ± 0.06 
7.46 ± 0.06 
5.59 ± 0.06 



Note. — All fluxes have been measured 
from our Gemini spectrum. No reddening eor- 
rection has been applied. There is an addi- 
tional systematic uncertainty in the flux mea- 
surements of ~ 10% due to fiux calibration. 

luminous. The peak radio luminosity observed for GRB 
100316D was a factor of a few higher than the first EVLA 
non-detection of SN2010ay (jMarguttil lin preparation! ) . 
The only relativistic explosion with detected radio emis- 
sion below o ur EVLA limits is the w eak and fast fading 
XRF 060218 (ISoderberg et al.l[2006al) . 

3. INITIAL CONSTRAINTS 

3.1. Light Curve Modeling 

We construct an R-band lightcurve for SN 2010ay us- 
ing the observations described in §2.31 (Table [2|). We 
convert the zpi and r' band data point s to the i?— band 
using t he fiducial lightcurve method of ISoderberg et ahl 
((20063) . assuming the unextincted i' — R and r' — R col- 
ors observed for the Type Ic-BL SN 1998bw at the ap- 
propriate epoch (|Galama et al.|[l998f) . Photometry for 
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Fig. 3. — Optical spectra of SN2010ay from Gemini (24 days after R-band peak) and the WHT (14 days). The spectrum is plotted 
decomposed into SN (above, with narrow lines clipped) and host galaxy (below, from Gemini, with splinc-fit subtracted) components for 
clarity The spectrum of SN 2010bh from LChornock et al.l 1I20T0I ') is given in black for comparison, at 21.2 da ys after the GRB 1 00316D 
trigger (~ 10.0 days after R-band peak, E^no et al. 201l3). The spectrum of SN 1998bw at +19 days from IPatat et al.1 l|20011 ) is also 
plotted. Both are transformed to the redshift of SN 2010ay. The SNe are shifted in flux for clarity. In the lower plot, relevant host galaxy 
emission lines are marked with a red line and labeled. 



the unfiltered CRTS images was reported bv lDrake et alj 
(|201Qi) after transformation to the synthetic V-band 
(A.J. Drake, private communication); we therefore con- 
vert to the i?— band assuming the V — R color of SN 
1998bw at the appropriate epoch. For the Pan-STARRSl 
photometry, the host galaxy flux was subtracted using a 
template image. For all other photometry, wc have sub- 
tracted the flux of the host galaxy numerically assuming 
the magnitude reported in Table [TJ A total (Galactic + 
host) reddening of E(B-V)= 0.2 mag has been assumed 
(see a2XT|l . 

To estimate the explosion date of SN2010ay we have 
fit an expanding fireb all model to the optical light curve 
(Figure [5]), following iConlev et al.l (|2006il . In this model, 
the luminosity increases as 



L oc 



t-tp 
1 + z 



(1) 



Wc derive an explosion date tp of 2010 February 
21.3±1.3. Here we have assumed an index n = 2. This 
suggests that the PSl 3n ipiband pre-discovery detec- 
tion image of SN2010ay was taken ~ 4 days after the 
explosion. Observations by the PSl survey have there- 
fore provided a valuable datapoint for estimating the ex- 
plosion date and also for constraining the rise-time of 
SN 2010ay, as well as other nearby SNe (e.g. SN2011bm 
iValenti et aTMy]3 ). 

In Figured! we compare the light c urve of SN 20 1 0ay to 
the SN Ib/c light-curve template of iDrout etldl (|2011[ ) 
after stretching by a factor of (1 -|- z). The template 
provides a reasonable fit to the optical evolution of 
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Fig. 4. — EVLA upper limits for SN2010ay (black arrows) arc compared with off-axis G RB afterglow lightcurve models (black curves; 
30, 60, 90°) and the observed radio light-curves for ordinary SNe Ib/c (grey: [S"od crbere"20071 and references within) and the radio afterglows 
of all GRB-SNe w ithin z < 0.25. SN2010ay is a fact or of 10^ to 10^ les s luminous than XRF 20903 (orange; ISoderberg et al.|[2004al'l . 
GRB 030329 (blue:IBerger et al.l2003bl : IFrail et al .'2005). and GRB 031203 llSoderberg et al.l2004bl V Relativistic, engine-driven SNe 1998bw 
Cred: IKulkarni et a .111998 1 and 2009bb (dark blue; Sodcrberg et al. 20101 ') are a factor of 10 more luminous than the SN 2010ay limits on a 
comparable timescale, while XRF 060218 lies a factor of a few below the limits. We constrain the radio counterpart to be no more luminous 
than XRF 060218 and comparable to the peak luminosities of ordinary SNe Ib/c. 



SN2010ay. Fitting the template to our photometry we 
derive (reduced = 0.3) a date of R-band peak of 2010 
March 18 ± 2 UT (MJD 55273 ± 2) and a R-band peak 
magnitude of M/j w —19.7 mag before extinction cor- 
rection. As discussed in §2.3. H based on the Balmer 
decrement observed for the host galaxy emission lines, 
we assume an extinction of E{B — V) = 0.2 mag. Ap- 
plying this correction, the peak absolute magnitude is 
Mr « -20.2 ± 0.2 mag. We note that this fitted value 
is « 0.2 mag fainter than that estimated from the data 
point near peak. Here, the uncertainty is dominated by 
the template fitting. 

Regardless of the extinction correction, SN 2010ay is 
more luminous th a n all the 25 SNe Ib/c in the sam- 
ple of IDrout et all (|20T1 . except for SN 2007D {Mr « 
—20.65 mag, which was also significantly extinctcd; 
Ay ~ 1.0 mag). Assum ing an intrinsic V-R color of zero 
at pe ak (e.g. 1998bw: ICalama et all [19981 iPatat et al.l 
[20?m . SN2010ay is also more luminous than any of 
the 22 GRB and XRF-pr oducing SNe in the compila- 
tion of iCano et al.l (|2011b[ ). all corrected for extinction, 
and is 1.5 standard deviations from the mean luminos- 
ity. The peak magnitude is only ^ 1 mag below that 



of the Type Ic SN2007bi (Mb = -21.3 ± 0.1 mag), 
which [Gal- Yam et al.l ()2009() report as a candidate pair- 
instability supernova. 

3.2. Large Nickel Mass for SN2010ay 

We use the available photometry for SN2010ay dis- 
cussed above to derive the mass of '"'^Ni required to power 
the optical light-curve under the assumption that the 
emission is powered by radioactive decay. Using the 
relatio n between AfNi and Mr^ found by IDrout et al.l 
(|20T1 . log(MNi) = (-0.41 Mr - 8.3) Mq, we esti- 
mate that SN 2010ay synthesized a nickel mass of Mni = 
0.9lg i[ Mq. Wc have estimated the uncertainty in the 
Mm estimate by propagation of the uncertainty in the 
template fitted peak magnitude — systematic uncer- 
tanties are not included. If we instead adopt the most 
luminous individual data point in the light curve as the 
peak value, rather than the smaller peak value from tem- 
plate fitting, we find Mm ~ 1.2 Mq. 

The Mm estimate for SN2010ay is larger than that 
of all but o i ie (SN 2007D) of the 25 SNe Ib/c in the 
IDrout et al.l ()2011f ) compilation and significantly larger 
than the estimate for GRB-SN 2010bh, Mni = 0.12 ± 
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O.OIMq (jCano et al.ll2011aD . On the other hand, Mni 
of SN2010ay is at least 3 x smaher than for SN2007bi 
{Mm « 3.5 - 4.5 Mq, lYoung et aLI [2010l) . A pair 
instabiUty supernova should produce Mni > 3 Mq 
(IGal-Yam et al.l[2009l) . 

3.3. Unusually high velocity 

As illustrated in Figure [3l the broad, highly-blended 
spectral features of SN2010ay at the time of the WHT 
observations (14 days after R-band peak, see §3.1|) re- 
sembl e those of SN 2010bh a t a similar time (10.0 d after 
peak. [Chornock et al.ll2010() . In particular, the blueshift 
of the feature near 6355 A is larger than in SN 1998bw 
and more similar to SN2010bh. This feat ure is com- 
monl y associated with Si II A6355A (e.g. iPatat et"al] 
|2001| ). However, this feature has two clearly-detectable 
absorption minima in SN 2010ay, but not in SN 2010bh. 
This could be due to increased blending in SN2010bh 
or the absence of contaminating lines. The red ends of 
the SN2010ay and SN2010bh spectra (rest wavelength 
> 7500 A) have similar P-Cygni features, but the emis- 
sion and absorption components in SN 2010bh are each 
blueshifted by ~ 200 A mo re th an in the spectrum of 
SN2010ay. iChornock et al.l ()2010f ) attribute this feature 
to the Ca II NIR triplet, with a (^/-weighted line cen- 
troid of 8479 A, and find a velocity that is high, but 
consistent with the early-time velocity of Si II A6355A 
(30 - 35 X 10^ kms-i). 

We measure the absorption velocity from the minimum 
of the Si II A6355A absorption feature (usi)- We smooth 
the spectr um using an inverse -variance-weighted Gaus- 
sian filter (iBlondin et al.l 120061 with dX/X = 0.005) and 
measure the minimum position of the redmost compo- 
nent of the absorption profile. The blue component of 
the absorption profile shifts blueward over time, suggest- 
ing that it is produced by a combination of ions separated 
by several 10'^ km s~^, such as He I A5876A and Na I D, 
whose relative optical depth changes with time. 

The absorption velocity inferred from the Si II A6355A 
feature is ~ 2x faster than that of SN 1998bw at simi- 
lar times, and more simil ar to that of SN 20 10bh (Fig- 
ure 3). For SN 1998bw, IPatat et al.l (|2001[ ) measured 
~ 10 X lO'^ km s~^ at +13 days. For SN2010bh, 
IChornock et al.l (|2010f ) measuredj^i w_26x 10^ km s^^ at 
-1-10.0 days after explosion. iPrietoi (|2010[ ) reported a ve- 
locity ofusi «22.6xl03 kms-i from the Si II A635 5A fea- 
ture in a spectrum of SN 2010ay taken at +4 days; iPrietol 
(|2010t) do not discuss the details of their methodology 
for the velocity measurement. From our [WHT,Gemeni] 
spectra taken [+14, +24] days after R-band peak (see 
^5+|) . we estimate vsi «[19.2,18.3] x 10^ km s'^. 

In addition to the broadening of the spectral fea- 
tures and the blueshift of the Si II A6355A line, ad- 
ditional lines of evidence suggest a high photospheric 
expansion velocity for SN2010ay. We measure vsi ~ 
[21.7, 20.1] X 10^ km s'^ from the Ca II NIR triplet on the 
smoothed [WHT, Gemini] spectra, relative to a line cen- 
ter at 8479 A. This is within a few 10^ km s~-^ of the usi 
we measure from Si II A6355A at these epochs. Further- 
more, we do not detect the broad emission "bump" near 
4500 A in either of our spectra. This feature was also 
absent in SN2010bh, but was identified in SN2003dh, 



TABLE 4 
Velocity Evolution of SNe Ic-BL 



SN 



SNe Ic-BL 



1997of 


6±2 


-0.8 ±0.4 


2002ap 


5±2 


-1.5 ±0.7 


2003jd 


10 ±1 


-0.5 ±0.2 


2007bg 


9±2 


-0.3 ±0.2 


2007ru 


10.3 ±0.7 


-0.5 ±0.1 


2010ay 


21 ±2 


-0.4 ±0.3 


Engine-driven SNe Ic-BL 


1998bw 


10.3 ±0.7 


-0.86 ±0.08 


2003dh 


12.0 ±0.9 


-0.9 ±0.1 


20031w 


10 ±1 


-0.8 ±0.2 


2006aj 


15 ±3 


-0.3 ±0.2 


2009bb 


10.7 ±0.4 


-0.84 ±0.08 


2010bh 


24 ±3 


-0.22 ±0.07 



Note. — To the velocity measurements for each SN, 
we have fit a power law of the form VQi — i;|^(t/30)", 
where t is the time since explosion in days and Vg^ is the 
velocity at 30 days in units of 10^ km s~^. These power 
law fits arc illustrated in Figure [5] Wc have assumed 
10% uncertainties in all velocity measurements. Uncer- 
tainties in the explosion date for SNe without detected 
GRBs vary due to availability of early-time photome- 
try; wc ad opt the followin g conservativ e uncertainties: 
7 d ( 1 99 7ef , |Hu et al.!|1997l) , 7 d (2002a p , IGal-Yam eTaP 

f200l^ . 6 d (2003jd, Valcnti craUl200gr . 14 d (2007bg) 
3d (2 007ru, |Sahu ct al..200al . 1 d (2009bb. [Pignata et"al1 

\201% , 2 d (2010ay; this paper) . 

SN2006aj, and seve ral Ic-BLs not as s ociate d with GRBs 
and normal SNe Ic; IChornock et al.l (|2010[ ) suggest that 
the absence of this feature indicates a high expansion ve- 
locity if it is due to blending of the iron lines to the blue 
and red of 4500 A. 

We compare the late-time expansion velocity of 
SN2010ay to other SNe Ic-BL and GRB-SNe with de- 
tailed, multi-epoch velocity measurements from the lit- 
erature in Figure [SJ We note that the velocities for SNe 
1997ef, 20 03dh, 20031w plotted in the figure were esti- 
mated by iMazzali et al.l ()2006D via spectral modeling, 
rather than measured directly from the minimum of the 
Si II A6355A absorption feature; however, these veloci- 
ties should be equivalent to within a few 10"^ km s~^, as 
the minimum of the Si II feature is typically well-fit by 
the photospheric vefocity of the sp ectral models (see e.g. 
iMazzah et alll2000l : iKinugasa et al., ,2002, '). In this figure, 
we also fit power-law gradients to the time-evolution of 
the velocity of these SNe with the form vsi = Wg?(i/30)", 
where t is the time since explosion in days and Wg? is the 
velocity at 30 days in units of lO'^ km s~^. These pa- 
rameters are listed in Table [H Figure [5] illustrates that 
most SNe are well described by a single power lawF^ 
However, due to lack of late-time spectroscopy, the Vgf 
measurement amounts to an extrapolation for some ob- 
jects (particularly 2006aj), and contamination from dif- 
ferent ions or detached features will add uncertainty to 
velocities measured from the Si II A6355A feature. 

SN 2010ay and 2010bh share high characteristic veloci- 

22 A break at ~ 2 X 10"^ km appear s to ex ist for SN 1998bw 
at ~ 16 days, as noted by IKinugasa et al] l|2002l 'l. 
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ties at 30 days after explosion and velocity gradients that 
are low relative to other broad-lined Ic SNe with and 
without associated GRBs. For SN 2010ay, u|P = 21 ± 2 
is 2 — 4 X larger than for other SNe Ic-BL without associ- 
ated GRBs = 6±2 for 1997cf, 5±2 for 2002ap, 10±1 
for 2003jd, 9 ± 2 for 2007bg, and 10.3 ± 0.7 for 2007ru) 
and is similar to the GRB-SN 2010bh (t>|0 = 24 ± 3). 



No other GRB-SN or SN Ic-BL has 



> 15. The 



SNe Ic-BL and GRB-SNc with the most shallow veloc- 
ity gradients among these twelve objects have a < —0.5; 
they are SN2006aj (a = -0.3 ± 0.2), SN2007bg {a = 
-0.3 ±0.2), SN 2010ay {a = -0.4 ±0.3), and SN2010bh 
{a = -0.22 ± 0.07). This places SN2010ay in the 
company of two GRB-SNe in having a slowly-evolving 
absorption velocity, and SN2007bg (whose unusually- 
fast decline rate d istinguishes it from other SNe Ic-BL, 
lYoung et al1l2010( ). The velocities of SNe 2010ay and 
2010bh, respectively, decline about 2 and 4x more slowly 
than the other SNe Ic-BL (mean and standard devia- 
tion: a = —0.8 ± 0.4) and about 1.5 and 3x more 
slowly the other GRB-SNe (a = -0.7 ± 0.2). The slow 
Si 11 A6355A absorption velocity evolution of SN 2010ay 
at late times resembles the slow evolution of the Fe II 
lines of the spectroscopical ly-normal Typ e Ic SNe 2007gr 
and 201 Ibm at late times (jValenti et al.ll2012[ ). 

Given the high peak luminosity of SN 2010ay fi )3.ip . we 
also consider the velocity of the candidate pair- instability 
SN 2007bi. Velocity measureme nts for SN 2007bi ar e only 
available at late times (> 50 d. lYoung et al.|[2010D . Fit- 
ting to these late-time Si II A6355A velocity measure- 
ments, we find that SN 2007bi has a characteristic ve- 
locity ~ 3x smaller than 2010ay: = 8 ± 2 and 
the late-time velocity gradient is ~ 2x more shallow: 
a = -0.2 ±0.2. 

3.4. Ejecta Mass and Energy 
We use the scaling relations provided bv IDrout et al.l 



()2011|) , based on the original formali sm of lArnettT i 1982[ ) 
and modified bv lValenti et al.l (|2008[ ). to derive the total 
mass of the ejecta and the kinetic energy. 



We assume the fitted peak magnitude for SN2010ay 
(Mr = -20.2 ± 0.2, ^31^) . the absorption velocity we 
measure from our WHT spectrum at 14 days after maxi- 
mum light {vsi — 19.2 x 10^ km s~-^, §3.3p . and a charac- 
teristic time (light-curve widt h) consistent with the data 
and the mean value from the IDrout et al.l (|201lD sample 
of SNe Ic-BL (r^ = 14 d). 

Using these values, the total mass ejected was Afoj ~ 
4.7 Mq, and the total kinetic energy of the explosion 
was Ek ~ 10.8 X 10^^ ergs. Hereafter we refer to the 
definition Ek^5i = Ej^/ 10^^ ergs. 

The systematic uncertainties associated with this mod- 
eling dominate the statistical uncertainties. In particu- 
lar, the models rely on the assumptions of homologous 
expansion, spherical symmetry, all ^^Ni centralized at the 



center of the ejecta, optically thick ejecta and constant 
opacity. 

We note that an earlier measurement of the absorption 
velocity is preferred for optical modeling. Since we have 
argued that SN 2010ay and SN 2010bh have similar char- 
acteristic velocities (wg-*), if we instead adopt a higher 
vel ocity of 25.000 krn s~^ as measured for SN2010bh 
bv ICano et al.l (I2011af ). we estimate Moj « 6.1 Mq and 
Ek,51 ~ 23.9 for SN2010ay. 

The Mcj and Ek. 51 of SN 20 10ay are c onsistent with 
the mean for SNe Ic-BL in the IDrout etlil.. ( ,2011 ) sam- 



ple {4:.7_^'g Mq and 11_4, respectively), because the au- 
thors assumed a velocity {vsi = 2 x 10''' km s""'^) similar 
to the late-time velocity we measure. For comparison, 
SN 2010bh had a total ejecta mas s of ~ 2 Mq and a 
total kinetic energy of Ek,5i ~ 13 (jCano et al.ll2011a[ ). 

The ratio of Ni to total ejecta mass is ~ 0.2 for SN 
2010ay. significantly higher than the values typical of 
SNe Ic-BL and GRB-SNc. For comparison, the ratio is 
just ~ 0.05 for SN 2010bh. Adopting the va lues derived 
from b olometric light curve modeling by iCano et al.l 
(|2011a[) . the Mm and Mm/M^j ratios for other GRB-SNe 
are: - 0.5 Mq and - 0.06 - 0.22 (1998bw), - 0.4 Mq 
and - 0.08 (2003dh), - 0.15 Mq and - 0.07 - 0.1 
(2006aj), and - 0.2 Mq and - 0.06 (2009bb). This 
ratio for S N 2010ay is l a rger t han that of all but 4 of the 
25 SNe of IDrout et all (|20T1 : the Type Ic SNe 2004ge 
(MNi/Mcj ~ 0.4) ans 2005eo {Mm/M^j ~ 0.2), the Type 
lb SN2005hg (- 0.4), and the Type Ic-BL SN2007D 
(^ 0.6). 

The large value of A/ni wc estimate for SN 2010ay raises 
the question of whether a process other than Ni decay 
may be powering its light-curve fe.g. lChatzopoulos et al.) 
I2009f ). An independent test of the physical process 
powering the light-curve is the decay rate of the late- 
time light curve which should be w 0.01 mag day~^ 
for SNe pow ered by ra dioactive decay of ^^Co. For 
SN 2007bi. iGal-Yam et al.. (,2009t ) derive Mm = 3.5 Mq 
from the measured peak magnitude and find that the 
late-time light curve is consistent with the decay rate 
of 56Co. While the Pan-STARRSl Stt survey also ob- 
served the field in 2011 March, the SN was not detected 
in our subtracted images and the limits are not constrain- 
ing in the context of ^®Co decay (sec H2.2I and Table 
[2). Another possible process is a radiation-dominated 
shock that emerges due to interaction with an opaque 
cir cumstellar mediuin, as h as recently been proposed 
by IChevalier fc Irwinl (120111) for the ultra- luminous SNe 
2006gv and 2010gx (jPastorello et al.ll2010l : IQuimbv et all 
1201 li) . However, while this class of ultra-luminous ob- 
j ects shares sonie spe ctroscopic similarities to SNe Ic 
(jPastorello et al.l [20100 . they sho w peak luminosities ^ 
4 - 100 X higher th an SN2010ay (|Chomiuk et all 120111 : 
IQuimbv et al][20Tl . 

4. CONSTRAINTS ON RELATIVISTIC EJECTA 

We use our EVLA upper limits for SN2010ay span- 
ning At w 29 — 433 days to constrain the properties 
of the Shockwave and those of the local circumstellar 
environment. The radio emission from SNe Ib/c and 
GRBs is produced by the dynamical interact ion of the 
fastes t ejecta with the surrounding material (|Chevalieii 
119820 . The kinetic energy of the ejecta is converted, in 
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Fig. 5. — A comparison of the time-evolving absorption velocity of SN2010ay and other SNe Ic-BL (red) and engine-driven explosions 
(blue) from the literature. For each SN, we fit a power law of the form tigi = Vg^{t/30)°', where t is the time in days, Ug? is the velocity 
at 30 days after explosion (dashed vertical line), and a is the velocity gradient. The velocities for SNe 1997ef, 2003dh, 20031w are from 
IMazzali et al.l I I2006I) . as determined by spectral modeling. The velocities for all other S Ne are measured from the Si II A6355A feature as 
follows: SN 2007ru are from[Sahu ct al. (20M); SNe 1998bw, 2006aj, and 2010bh are from lChornock et al.l |[20TOf), fr om spectra in references 
therein; SN 2007b g are fr om Young ct al. (2010); SNe 2002ap, 2009bb, and 2003jd are from IPignata et al.l idOUT) and references therein; 
SN 2010ay is from lPrietcl | |2010I) and this paper. See Table l4l for uncertainties. 



part, to internal energy of the shocked material which it- 
self is partitioned between relativistic electrons (ce) and 
amplified magnetic fields (es). Following the breakout 
of the Shockwave from the stellar surface, electrons in 
the environment of the explosion are shock-accelerated 
to relativistic velocities with Lorcntz factor, 7e and dis- 
tributed in a power-law distribution characterized by 
Ni^je) oc 7~^. Here, p characterizes the electron en- 
ergy index. The particles gyrate in amplified magnetic 
fields and give rise to non-thermal synchrotron emission 
that peaks in the radio and mm-bands in the days to 
weeks following explosion with observed spectral index, 
Fi, (X At lower frequencies the emission is sup- 

pressed due to sync hrotron self-abso rption which defines 
a spectral peak, Vp (jChevalierl 1 1 998f ) . 

The dynamics of the Shockwave determine the evo- 
lution of the synchrotron spectrum, and in turn, the 
properties of the observed radio light-curves. In the 
case of SNe Ib/c, there are three primary scenarios for 
the dynamical regime of the ejecta depending on the 
shock velocity, v = Pc (associated Lorentz factor, F): 
(i) non-relativistic {v ~ 0.2 c) free-expansio n as in the 
case of ordinary SNe Ib/c (jChevalieii |1998() . (ii) a de- 
coupled and relativistic (F ~ 10) shell of ejecta that 
evolves accordi ng to the Blandf ord-McKee solution for 
several months (jSari et al.lll998| ) before transitioning to 



the Sedov- Taylor regime (jFrail et al.ll2000( ). This is the 
standard scenario for typical GRBs. And (iii) a sub- 
energetic GRB with trans-relativistic velocity (/3F < 3) 
that bridges the free-expansion a nd Blandford McKee 
dynamical regimes (e .g. SN 1998bw; iKulkarni et al.l[T998l 
ILi fc Cheval"ieH[l999l) . 

We consider our EVLA upper limits in the context 
of these three models below. For shock velocities of 
y > 0.2c, Ee ~ 0.1 is reaso nable (jSoderberg et al.l [20051 : 
iChevalier fc Franssonll2006l ) . Wc further assume equipar- 
tition, Ce = cb = 0.1. We adopt a free expansion model 
for both the non-relativistic ordinary SN Ib/c case and 
the sub-e nergetic, trans- r elativ istic GRB scenario. As 
shown by ILi &: Chevalieii ()1999[ ) a free-expansion model 
is still reasonable in the trans-relativistic regime (cases 
iii, see above). 

4.1. Freely- expanding shockwave 

In the free-expansion scenario, a shock discontinuity 
separates the forward and reverse shocks, located at the 
outer edge of the stellar envelope. The bulk ejecta is in 
free expansion while the thin layer of post shock material 
is sli ghtly decelerated, R cx t^'^ () Chevalier fc FranssonI 
[20061) . At a given frequency, the bell-shaped light-curves 
of the SN sync hrotron emission may be described as 
(|ChevaliedH"998[) 
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1.582 X L^^p 



At 



1 



(At/tp) 



~(o + i>) 



(4) 



where L^.p is the flux density at the spectral peak 
at epoch, tp. Assuming an electron index of 
p K, 3, consistent with r adio spectra of SNc Ib/c 
(|Chevalier fc Franss on '2006*). the exponents are a w 2.3 
and h « 1.3. The time averaged Shockwave velocity is 
V « R/ At where R is the Shockwave radius defined as 



i?«2.9 X 10 



16 



-1/19 



1028 erg s-i Hz- 



9/19 



V5 GHz/ 



(5) 



Here we make the assumption that the radio emitting 
region is half of the total volume enclosed by a spheri- 
cal blastwave. Next, we estimate the internal energy, E, 
of the radio emitting material from the post-shock mag- 
netic energy density, E « B^R^ /\2eB where w e maintain 
the as sumption of equipartition. As shown by iChevalierl 
()1998[) . the amplified magnetic field at peak luminosity 
is also directly determined from the observed radio prop- 
erties, 



B«0.43 
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.5 GHz; 

Finally, the mass loss rate of the progenitor star, 
M , may be derived from the number density of emit- 
ting electrons. Here we normalize the wind profile ac- 
cording to p cx Ar~2 and A* = A/5 x 10^^ g cm~. 
This normalization of A* implies that an A* of 1 corre- 
sponds to typical Wolf-Rayet progenitor wind properties 
of M ~ 10^^ Mq yr^^ and a progenitor wind velocity of 
u,„ = 10^ km s^^ 



A, wO.15 



VO.l/ 





-8/19 ^ 









1028 erg s-i Hz-i 
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At 



5 GHz/ V 10 days 

where we assume a shock compression factor of ^ 4 and 
an nucleon-to-electron ratio of two. 

We built a two-dimensional grid of fiducial radio light- 
curves according to Eqn. |4] in which we vary the pa- 
rameters Li/p and Vp over a reasonable range of param- 
eter space, bounded by tp « [1,3000] days and F^^p « 
0.04 — 1000 mJy. We identify the fiducial light-curves as- 
sociated with a radio luminosity higher than the EVLA 
upper limits for SN2010ay at each epoch as these are 
excluded by our observations. We extract the physical 
parameters associated with these excluded light-curves 
(i?, B, E, A^) to define the parameter space excluded by 
our radio observations. The parameter space for Vp and 
F^^p arc bounded by the respective values for which the 
model exceeds relativistic velocities, /3r ~ few. 



(7) 



As shown in Figure [SI our deep EVLA limits enable us 
to rule out a scenario in which there is copious energy, 
E > 10*8 gj.g^ coupled to a relativistic outflow, in this 
two-dimensional E — v parameter space. The excluded 
region includes GRB-SNe 1998bw and 060218 as weU as 
the relativistic SN 2009bb. It docs not exclude the stan- 
dard scenario in which a small percent of the energy is 
coupled to fast moving material within the homologous 
outflow, as is typically observed for ordinary SNe Ib/c 
{E « lO''^ and V « 0.2c; ISoderberg et"all[2010l) . 

Next we consider the effects of circumstellar density 
since lower mass loss rates produce fainter radio coun- 
terparts. As shown in Figure [71 the EVLA limits for 
SN2010ay exclude the region of parameter space popu- 
lated by SNe 1998bw and 2009bb with mass loss rates 
of A* ~ 0.1, however, the low density environment of 
GRB 060218 lies outside of our excluded region due to 
its lower GSM density, ^ 0.01 that gives rise to a 
lower luminosity radio counterpart. 

4.2. Relativistic Ejecta 

In the case of relativistic deceleration the ejecta are 
confined to a thin jet and are physically separated from 
the homologous SN component. Deceleration of the jet 
occurs on a timescale of At » (Esl /A^,) years in a 
wind-stratified medium (iWaxmanI l2004f ) . On this same 
timescale, any ejecta components that were originally 
off-axis spread sideways into the observer's line-of-sight. 
While the early EVLA limits constrain the properties of 
the on-axis ejecta according to the free-expansion model 
described above, the late time EVLA upper limits con- 
strain any radio emission from a GRB jet originally 
pointed away from our line-of-sight. 

For this scenar io, we adopt the semi-analytic model of 
ISoderberg et all (pOedi ) for off-axis GRB jets. In Fig- 
ure [4| we compare the radio upper limits for SN 2010ay 
with the predictions for an off axis jet with standard 
parameters {Ek ~ 10^^ erg, A, = 1, = es = 0.1, 
6j = 5°) and assuming a viewing angle of 9oa = 30, 60, or 
90 degrees. Our EVLA upper limits rule out all three of 
these model light-curves. We derive the two-dimensional 
parameter space (energy and CSM density) that is ex- 
cluded based on our EVLA upper limit at A< « 1.2 years. 
We note that this model accommodates the full transi- 
tion from relativistic to non-relativistic evolution. We 
built a collection of model light-curves spanning param- 
eter range, A* « [0.01 - 100] and E w [10^^^ - 10^2] ^j-g^ 
maintaining the assumption of p = 3 and = es = 0.1. 
Here we adopted a jet opening angle of 9j = 5° and an 
off-axis viewing angle of 6'oa = 90° (the most conserva- 
tive s cenario) . We note that Ivan Eerten fc MacFadvenI 
(|2011| ) have developed off-axis GRB afterglow lightcurve 
models based on hydrodynamic simulations that repro- 
duce the semi-analytic n iodels presented here t o within 
factors of a few (see also Ivan Eerten et al1[20Tot) . 

As shown in Figure [51 we are able to exclude the pa- 
rameter space associated with typical GRBs, i.e. E w 
Iq5i gj.g (beaming corrected) and A* w 1. GRBs with 
lower energies and densities are better constrained using 
the trans-relativistic formalism above. In conclusion, our 
radio follow-up of SN 2010ay reveals no evidence for a rel- 
ativistic outflow similar to those observed in conjunction 
with the nearest GRB-SNe, however a weak afterglow 
like that seen from XRF 060218 cannot be excluded. 
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Fig. 6. — The region of energy- velocity space ruled out (red) by 
our EVLA observations for on-axis cjccta under the assumption of 
a free-expansion model. 
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Fig. 7. — The region of energy-mass loss space ruled out (red) by 
our EVLA observations for on-axis ejecta under the assumption of 
a free-expansion model. 

5. CONSTRAINTS ON AN ASSOCIATED GRB 

Given the estimate of the explosion date we have de- 
rived ( ^3.ip . we have searched for gamma-ray emission 
that may have been detected by satelhtes. No GRBs 
consistent with SN 2010ay were reported in the circulars 
of the Gamma-ray Coordinates Network Circulars, but 
it is possible that bursts were detected below the instru- 
ment triggering thresholds. 

We next consulted the sub-thres hold bursts from the 
Swift Burst Alert Tele scope (BAT; iGehrels et all 120041 : 
iBarthelmv et al.l |2005() detected within the 6 days sur- 
rounding the explosion date estimate. We find that 
no gamma-ray emission was detected within 0.5 degrees 
of the position of SN2010ay by the BAT during this 
period. Given the sensitivity of the BAT, this cor- 
responds to an upper limit on the gamma-ray flux of 
~ 10~^ erg s~^ cm~^ (15-150 keV). However, the field 
of SN2010ay was in the field of view of the instrument 
for only 106 ksec during these 6 days, or ^ 20% of the 
duration. 




Mass Loss Rale (A.) 

Fig. 8. — Our EVLA observations at At ?s 1.2 years after explo- 
sion constrain the properties of a possibly associat ed of f-axis GRB 
jet. Using our semi-analytic model as described in ^4.21 we assume 
partition fractions of ee = eb = 0.1, 9j = 5 degrees, p = 2.5, 
and an off-axis viewing angle of 6oa = 90 degrees. We are able to 
exclude the region of Ef^ — A* parameter space (red) associated 
with typical GRBs, i.e. Ejf = 10^^ erg and At = 1 (dashed black 
lines). 

For complete temporal coverage, we have searched the 
records of the interplanetary network (IPN), which is sen- 
sitive to bursts with fluences down to ~ 6 x 10"'' erg cm~ ^ 
(25-150 keV) (50% eSiciency hmit, IHurlev et alTl20Tot) . 
and observes the entire sky 

with a temporal duty cycle close to 100%. An unde- 
tected, 

sub-threshold burst should have a fluence below this 
limit. Between 2010 February 21 and 25, inclusive, a 
total of 12 bursts were detected by the spacecraft of the 
IPN (Mars Odyssey, Konus-Wind, RHESSI, INTEGRAL 
(SPI-ACS), Swift-BAT, Suzaku, AGILE, MESSENGER, 
and Fermi (GBM)). Ten of them are confirmed bursts; 
they were observed by more than one instrument on 
one or more spacecraft, and could be localized. Two 
of them are unconfirmed bursts; they were observed by 
one instrument on one spacecraft (Suzaku). The total 
area of the localizations of the 10 confirmed bursts is 
^ 0.58 X Att sr. This implies that about 0.58 bursts can 
be expected to have positions that are consistent with 
any given point on the sky simply by chance (i.e. within 
the 3(T error region) , and indeed none of the bursts in this 
sample has a position consistent with the SN position. 

These non-detections imply upper-limits to the 
gamma-ray energy (E^) of a burst that may have been 
associated with SN 2010ay. The IPN non-detection in- 
dicates <6x lO'*^ erg (25-150 keV), while the BAT 
non-detection indicates that the peak energy of the burst 
was < 1 X 10'^'' erg s~^ (15-150 keV) if the burst occurred 
while in the field of view of the instrument. 

We consider whether or not a hypothetical GRB as- 
sociated with SN 2010ay would have been detected by 
BAT or the IPN if it had characteristics similar to well- 
studied SN-associated GRBs. The isotropic prompt en- 
ergy release of long GRBs is typically E^^iso ~ 10^^ erg, 
however the prompt emission of the sub-energetic class 
of GRB-SNe can be several orders of magnitude fainter 
(|Soderberg et al.„2006M . GRB 980425/SN 1998bw had 
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a peak ga mma-ray luminosity of 5 x 10"*^ erg s ^ (24- 
1820 keV iGalamaeFaLl [19981 ) ■ which is a factor of 2 be- 
low o ur BAT hmit, and E-yAso ~ 5 x 10^'' erg (jPian et alJ 
|200(I ). more than a factor of five below our limit. Neither 
satellite should have detected such a burst. In contrast, 
GRB 031203/SN20031W had a peak gamma-ray lumi- 
nosity of 1 X 10"^ erg s-i (20-200 keV) and a total 
isotropic eq uivalent energy of E^ iso = (4± 1) x lO**^ erg 
(jSazonov et al. 200j), about two orders of magnitude 
above the sensitivity of the BAT and twice the threshold 
of the IPN, respectively. GRB 030329/200 3dh was even 
more luminous, with iso 7 x 10^^ erg (|Hiorth et al.l 
|2003[) . A burst like GRB 031203 or 030329 should cer- 
tainly have been detected by IPN, or the BAT if it oc- 
curred while the field of SN 2010ay was in the field of view 
of the instrument. XRF 060218/SN2006aj, an extremely 
long-duration {At « 2000 s) event, had a peak luminosity 
observed by BAT of - 10"^ erg cm^^ g-i (15.150 keV), 
corresponding to a peak emission of w 2.4 x 10^^ erg s~^ 
given the redshift of the burst (z = 0.033), and a to- 
tal isotropic equivalent energy of E'^^iso = (6.2 ± 0.3) x 
1q49 gj.g ([Campana et al.|[2006D . If such a burst was asso- 
ciated with SN2010ay, its peak emission may have been 
a factor of four below the BAT sensitivity limit, although 
its total isotropic energy emission is an order of magni- 
tude larger than our IPN limit for SN 2010ay. Finally, the 
event whose host galaxy and supernova properties seem 
most similar to SN 2010ay, GRB 100316D/SN 2010bh 
had E^^iso > (5.9 ± 0.5) x ere lStarling eFall ([20T1 
— a full order of magnitude above our IPN limit. 

Another possibility is that prompt emission associ- 
ated with SN2010ay may have been too soft to be de- 
tected by the BAT or IPN. For example, the spectrum 
of XRF 060218 rose to a peak at 0.3 - 10 keV at - 985 
s after triggering, then softened significantly thereafter. 
Even though the the total emission of this burst is well 
above our IPN limit, it may have escaped detection if it 
was similarly soft. 

6. SUB-SOLAR HOST ENVIRONMENT METALLICITY 

We estimate the oxygen abundance of the host environ- 
ment of SN 2010ay from the strong nebular emission line 
fluxes measured from our Gemini spectrum (Table[3|). At 
the distance of the host galaxy, the 1 " Gemini slit width 
corresponds to a physical size of 1.3 kpc. The properties 
we infer for the explosion site of SN2010ay represent a 
luminosity-weighted average over this radius. 

We employ several different oxygen abundance di- 
agnostics in order to determine the met allicity of the 
host galaxy from its optical spectrum fe. g.lMod iaz et al.l 
120 111) From the 03N2 diagnostic of IPettin i & Pagel 
(j2004[) (PP04), we derive a metalhcity of log(0/H) = 
8.19, or Z ~ 0.3 Z©, ad opting the solar me tallicity 
log(0/H)o-|-12 = 8.69 from lAsplund et al.l (l2005l) . Using 
the N2 diagnostic of PP04, we find log(0/H)-hl2 = 8.26. 
Using the abundance diagnostic, R23 = log([0 II]A3727-|- 
[O III] AA4959 , 5007) /H/3, we find log(0/H)+12z94 = 

8.49 (Za ritslwlFaLl [T99I a nd log(O/H)+12KD02 = 

8.50 pCewlev fc Dopital 120021 ). However, these R23- 
based estimates are more sensitive to flux-calibration 
and reddening-correction. Moreover, there is a well- 
known bi-valued relationship between R23 and oxygen 
abundance. The value R23 = 0.873 ± 0.006 measured 



at the explosion site places it near the turnover point, 
but we assume that it lies on the up per branch based 
on its [N II] /[O II] ratio, following iKewlev fc EllisonI 
(|2008| ). The metaUicity values we derive using the 
PP04 and KD02 diagnostics are approximately equiv- 
alent given the offset that e xists between these two 
diagnostics iKewlev fc EThso^ (gOO^. These measure- 
ments are similar to the values reported bv lModiaz et al.1 
(I20TI (log(0/H)-hl2 [PP04,KD02]= [8.2,8.4]) for the 
SN2010ay host galaxy. The statistical errors in our 
strong line metallicity estimates are small (< 0.01 dex), 
as determined by propagating the errors in the line flux 
measurement through the abundance calculation. How- 
ever, for example the representative systematic error for 
the PP04 03N 2 abundance d iagnostic is ~ 0.07 dex, as 
determined by IKewlev fc EU ison (200l) via comparison 
to other strong line abundance indicators. 

Fortunately, our detection of the weak [O III] A4363 
auroral line (S/N'--' 6; Figure 3) allows us to derive an 
oxygen abundance via the "direct," Te method. We em- 
ploy a^nethodology similar to that used by, for exam- 
ple, iLevesque et al.l (|2010a| ). We first derive the elec- 
tron temperature (Tg = 1.09 ±0.06 x 10** K) and density 
(rie = 80 ± 20 cm-3) from the [O III] and [S II] line ra- 
tios using the temden task of the IRAF package nebular 
(jShaw fc DufouHll994D. der i ve the 0+ temperature using 
the calibration of iGarnet^ (|1992[ ). an d finally e s timat e 
the 0+ and 0++ abundances following IShi et al.l (|2006l ). 
The direct abundance, log(0/H)+12 = 8.24 ± 0.08, is in 
good agreement with the PP04 03N2 value. The stated 
uncertainty reflects the propagation of the uncertainties 
for the line flux measurements. Indeed, the offset be- 
tween thes e two diagnostics should be very small at this 
metallicity [Kewlev fc ElUsonl ()2008l ). 

We estimate the star formatio n rate (SFR ) of th e host 
galaxy using the Ha relation of iKennicuttI (jl99^ . Af- 
ter correcting for host galaxy extinction, we measure the 
Ha luminosity from our Gemini spectrum (Table [3]) and 
estimate SFR= 1.1 Mq yr-^ 

6.1. Blue Compact Galaxy Host 

We compare the host galaxy of SN2010ay to the 
nearby galaxy population of the SDSS spectroscopic sur- 
vey. The physical properties of the host galaxy, SDSS 
J123527.19-t-270402.7, are estimated in the MPA/JHU 
catalo ePl . The total (photometric) galaxy stellar mass 
(M*) is given as 3.6ti'3 ^ Mq, the aperture-corrected 
SFR is l.OtpJ Mq yr~\ and the nuclear (fiber) oxy- 
gen abunda nce (O/Hp) is l o g(0/H )-|-12 =8.58t|^:^3 on 
the scale of iTremonti et al.l ()2004l ) (T04). The specific 
star formation rate (SSFR) of the host galaxy is then 
~ 2.8^q'^ Gyr~^. For consistency, we consider these val- 
ues of M*, the oxygen abundance, and the SFR for the 
host galaxy of SN 2010ay when comparing to other galax- 
ies in the MPA/JHU catalog. 

The oxygen abundance and SFR of the host galaxy 
of SN2010ay listed in the MPA/JHU catalog are con- 
siste nt with the values we derive in this paper (see 
also iKellv fc KirsWl l20Tl[) . The MPA/JHU cata- 

■^^ http: //www. mp a -garching.mpg.de/SDS S (described in 
|K.aufrmann ct al." '200% ITremonti etall l200l : IBrinchmann et al.l 
12004 : .SaMuQt al,, 200^, and updated for SDSS DR7) 
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log lists metallic i ties o n the T04 scale. Using the 
iKewlev fc EUisonl (I2008D conversion to the PP04 scale, 
the T04 metallicity estimate corresponds to a metallic- 
ity of log(0/H)+12 = 8.38, which is - 0.2 dex higher 
than the one we measure (log(0/H)+12 = 8.19). How- 
ever, there is a large 0.2 dex) rms scatter between 
the PP04 03N2 and T0 4 diagnostics at the reg ime of 
log(0/H)+12ppo4 ~ 8.2 (IKewlev fc EllisonI 120081) . The 
SFR in the MPA/JHU catalog is also in good agreement 
with the value we estimate from the Ha luminosity. Al- 
though our estimate does not include an aperture cor- 
rection, the size of the Gemini slit (1") should encom- 
pass most of the star formation in the galaxy (Petrosian 
r = 1.355"; Table [T]). 

The mass to light ratio of the host galaxy of SN 2010ay 
is low compared to typical star-forming galaxies. To 
compare the host galaxy to the general galaxy popu- 
lation, we select a subset of the MPA/JHU catalog by 
requiring that estimates of , SFR, a nd O/Hp be avail- 
able a nd we remove AGN according to iKauffmann et al] 
(|2003f ). We consider 167,837 starbursting galaxies fol- 
lowing these constraints. The host galaxy ranks in 
the [4th,38th,llth] percentile in [M*,SFR,0/Ho] among 
these galaxies. Among the selected galaxies with a stel- 
lar mass as low as the host galaxjo, the median and 
standard deviation of the B-bancQ absolute magnitude 
is -15.8 ± 1.3 mag. With Mb =-18.35 ± 0.05 mag, the 
host of SN 2010ay is more luminous than other galax- 
ies with a similar mass at the 2a level. The discrepancy 
is due to the blue color of the SN2010ay host galaxy, 
which indicates a stellar population that is very young 
and therefore has a low stellar mass to light ratio. Among 
the 1,184 galaxies in the MPA/JHU catalog that meet 
the constraints above and have a color similar to the 
host of SN2010ay (0.47 <u-r < 0.67, from SDSS fiber 
magnitudes), the host galaxy has typical properties, with 
[A/,,SFR,0/Ho] in the [46th,49th54th] percentile. 

Based on these properties, we classify the host galaxy 
of SN 2010ay as a luminous Blue Compact Galaxy 
(BCG). BCGs span a large range in luminosity (—21 < 
Mb < —12, where luminous BCGs have AIb < —17), but 
are distinguished by their blue colors {B — V < 0.45), 
high SFR (1 < SFR < 20 Mq yr'^), and low metal- 
licity (Z(^/50 < Z < Zq/2; iKunth fc OsTIh] [2OOOI : 
IKong fc CheMl[200l . The host galaxy of SN2010ay 
has a luminosity {Mb =-18.35 ± 0.05), color {B - V ^ 
0.11 ± 0.07), SFR {l.Otoi Mq yr'^), and metallicity 
{Z ~ 0.3 Zq) consistent with all these ranges 13 

6.2. Comparison to SNe Ic-BL and GRB-SNe Host 
Galaxies 

Our measurement of the metallicity from the Gemini 
spectrum indicates that the explosion site of SN 2010ay 
is ~ 0.5 (0.2) dex lower in metallicity than the median 

^"^ This subset is selected such that the host galaxy of SN 2010ay 
has the median mass: 1.39 < M* < 5.62 X 10* Mq, N^i^ = 6, 978. 

We obtain B-band magnitudes by converting the fc-corrected 
gri magnitudes given in the MPA/JHU catalog t o BVR magnitudes 
using the transformation of lBlanton fc RoweisI 1)20071 ) . 

We note that a large fraction of luminous BCGs show evidence 
for disturbed morph o logies or interaction with close companions 
IIGarland et al.1 120041 : IL6pez-Sanchez et al.| 120061 ). but we do not 
see evidence for a companion at the depth of SE)SS images of the 
host galaxy of SN 2010ay. 



SNe Ic (Ic-BL) in the sample of lModiaz et all ()201lD . In 
that sample, the median PP04 03N2 metallicity mea- 
sured at the explosion site of SNe Ic is log(0/H)-|-12 « 
8.7 and for Ic-BL is « 8.4 dex, for 12 and 13 objects, 
respectively. If instead the KD02 metallicity is used, the 
median of the sample is « 8.9 dex for SNe Ic (13 objects) 
and w 8.7 dex for Ic-BL (15 objects), so the abundance 
of the SN2010ay host galaxy is similarly low compared 
to the median. 

The metallicity of the environment of SN 2010ay 
is more similar to previously-studied nearby GRB-SN 
progenitors. A metallicity identical to our measure- 
ment was measured a t the explosion site of SN2010bh 
(iLevesaue et al.ll2011l) : logfO/H)-H2 = 8.2. In the sur- 
vev of iLevesaue et al.l ()2010al) . and adding the measure- 
ment for SN 2010bh, the GRB-SNe host galaxies have an 
average and standard deviation PP04 03N2 metallicity 
of log(0/H)-f 12 = 8.1 ± 0.1 on the PP04 scale, which is 
consistent with the SN2010ay env ironment. Among the 
17 LGRB host galaxies surveyed in lSavaglio et al.l (|2009f ) 
the average metallicity is somewhat lower, 1/6 Zq or 
log(0/H)-f 12 ~ 7.9, but these are at an average redshift 
of z ^ 0.5 that is much higher than SN 2010ay. 

This evidence suggests that the host galaxy of 
SN2010ay has chemical properties more consistent with 
LGRBs/GRB-SNe than SNe Ic-BL without associated 
GRBs; however, selection effects may mitigate this dis- 
crepancy. SNe found in targeted surveys of luminous 
galaxies have host galaxy properties biased towards 
higher metallici ties, due to the lum inosity-metallicity 
(L-Z) relation (jTremonti et al.ll2004[) . LGRBs are found 
in an untargeted manner through their gamma-ray emis- 
sion and therefore are not biased by this relation. 

SN2010ay joins a growing list of SNe Ic-BL that 
have been discovered in low metallicity host galaxies. 
Given the systematic uncertainty in strong line oxygen 
abundance diagnostics (~ 0.07 dex), we will consider 
host galaxies with metallicity Iog(0/H)ppo4 + 12 < 8.3 
(Z < 0.4^0) to be in the low-metallicity regime of 
SN 2010ay. Among the 15 SNe Ic-B L (9 discovered by un - 
targ eted searches) i n the surveys of lModjaz et al.l ()2008l ) 
and iModiaz et all (|2011| ). 4 were found in low metal- 
hcity environments: SN [2007eb,2007qw,2005kr,2006nx] 
at Iog(0/H)ppo4 + 12 =[8.26,8.19,8.24,8.24]. AU of 
these SNe w ere d iscovered by untargeted searches. 
lYoung et al.l ^2(M ) measure the metallicity of the 
host galaxy of the broad-lined Ic SN2007bg to be 
Iog(0/H)ppo4 + 12 = 8.18, ahhough this SN has 
lightcurve and spectral properties that distinguish 
it from normal SNe Ic-BL (i |3.3p . Furthermore, 
lArcavi et al] (|2010l ) find that SNe Ic-BL are more com- 
mon in dwarf {Mr > —18) host galaxies, which the au- 
thors attribute to a preference for lower metallicities. 

The star formation properties of the host galaxy of 
SN2010ay also resemble the host galaxies of LGRBs. 
If we consider those galaxies in the MPA/JHU catalog 
with masses similar to the host galaxy of SN 2010ay 
(as defined above), then the median SFR and O/Hq of 
these galaxies is 0.13 Mq yr~^ and log(0/Ho)-|-12 = 
8.36, respectively. The host galaxy of SN2010ay is in 
the [96th, 77th] percentile for [SFR,0/Ho] among these 
galaxies. This indicates that, while the host galaxy of 
SN 2010ay falls within la of the mass-metallicity (M — Z) 
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relation for star-forming galaxies, its SFR is extreme for 
its mass. The 39 LG RB host galaxies in the survey of 
iSavaglio et alj (|2OO90 are similarly low in mass and have 
high star formation rates, with an average stellar mass 
of ~ 10^ Mq and SSFR- 3.5 Gyr"!. 

The host galaxy of SN2010ay falls below the L - Z 
relation for nearby star-forming galaxies, as illustrated 
by Figure [9l We have transformed the fc-corrected 
nri magnitudes from the MPA/JHU catalog to B-band 
(|Blanton &: RoweisI [20071 ). At the luminosity of the host 
galaxy of SN2010ay, the median metallicity and stan- 
dard deviation of the SDSS galaxies on the T04 scale is 
log(0/H)-(-12 =8.93±0.17; the host galaxy of SN2010ay 
falls in the 3rd percentile. In other words, the host galaxy 
of SN201 0ay is a 2a outlier fro m th e L — Z relation . 
Similarly, iLevesQue etHI (|2010a[) and IHan et all (|2010f ) 
suggest that the host galaxies of LGRBs fall below the 
L — Z relation as defined by normal star- forming galaxies, 
B CGs, and the ho s t gala xies of Type Ic SNe. 

iMannucci et al.l (|2011[ ) have explained the offset of 
LGRB host galaxies from the M — Z relation as a 
preference for LGRBs to occur in host galaxies with 
high SFR, as characterized by the Fundam e ntal M etal- 
hcit y Relation (FMR) of 



TABLE 5 

Comparison between SN2010ay and GRB 
100316D/SN2010BH 



Mannucci et al.l (|20100 (see 



also iLara-Lopez et al.l 120101 ). Using the extended FMR 



for low mass galaxies from (jMannucci "eFalllmTI) . the 
host galaxy of SN 2010ay should have a metallicity of 
log(0/H)-Hl2 = 8.20 give n its stellar inass an d SFR. The 
FMR is calibrated to the iNagao et all (j2006[) metallicity 
scale, which is similar to that of PP04 at this metallic- 
ity. Given the intrinsic scatter in the extended FMR on 
the order of ^ 0.05 dex, this value is consistent with 
the PP04 value we m easure from our G e mini spectrum: 
log(0/H)+12 = 8.19. IKocevski fc (pOTTI) similarly 

explain the offset of LGRB host-galaxies from the M — Z 
relation as a SFR effect, but suggest that the long GRB 
host galaxies have even higher SFR than would be im- 
plied by the FMR. 

SN 2010ay is an example of a SN Ic-BL where the host 
galaxy is consistent with the M — Z relation for star- 
forming galaxies, but deviates from the L — Z relation 
due to its low stellar mass to light ratio (t j6.ip . Its 2a 
discrepancy from the L — Z relation would be hard to 
explain as a SFR rate effect alone because among galax- 
ies in the MPA/JHU catalog without AGN (as defined 
above) and with Mb within 0.1 mag of the host galaxy 
of SN 2010ay, the host galaxy has a SFR in the 26th per- 
centile (< la discrepancy). 

7. DISCUSSION 

SN 2010ay has all the hallmark features associated with 
GRB-SNe, and yet we find no evidence of a relativis- 
tic explosion to sensitive limits. We are able to place 
constraints on the energy, density, velocity, progenitor 
mass-loss rate, and gamma-ray flux of any GRB that 
may have been associated with it. In particular, we may 
rule out the association of a GRB that looks similar to 
any spectroscopically confirmed GRB-SN to date, except 
for XRF 060218. 

The low metallicity of the host environment of 
SN2010ay may be suitable for GRB jet forma- 
tion in the "collapsar" model, but our observations 
strongly constrain any re l ativist ic outflow (ij4] & [5]). 
In iMacFadvcn fc WooslevI (|1999|) . a high rate of ro- 



Property 


SN 2010bh 


SN 2010ay 




Host galaxy properties 


log(0/H)-(-12 ^ 


8.2 


8.19 


Rodshift (z) 


0.059 


0.06717 


Mr 


-18.5 


-18.94 


Explosion properties 


SN type 


Ic-BL 


Ic-BL 


t)|P (10=* km s-i)'' 


24 ±3 


21 ± 2 


Mr 


-18.60 ±0.08 


-20.2 ±0.2 


Mm (Mq) 


0.10 ±0.01 




Mej (Mq) 


1.93-2.24 


Z 4.7 


Ek,51 


12.0-13.9 


Z 10.8 


GRB energy (Eiso, erg) 


> 5.9 X 10"^ 


< 6 X 10*8 



Note. — The observe d properties of SN2010 bh and its 
host galaxy are given bv IChornock et ah| (120101) and light 
curve modelling was performed by lCano et all I2011al) . The 
properties of SN2010ay are derived in tiiis paper. 
^ The oxygen abundance derived from the PP04 03N2 
metallicity diagnostic, as discussed in Jj6] 

^ The absorption velocity at 30 days after explosion, as mea- 
sured from the Si ll A6355A feature in 

^ The lower limit of the total isotropic energy release esti- 
mated bv lStarling et all 120111) . 



tation in the core of the progenitor is required to 
power a relativistic jet. A low metallicity is pre- 
scribed to suppress the line driven winds that would 
deprive the core of angul ar momentum . Ap parently 
supporting this model, iStanek et al.l (|2006l ) found 
that the isotropic prompt energy release of the GRB- 
SNe decreases steeply with metallicity, and other 
surveys have found observational evidence for the 
preferential occu rrence of GRB-SNe i n low-metallicity 
host 



20041 SoUcrman et al. 120051 Modiaz et al.l 


2006; 


Wiersema et al.l 120071 


Christensen et al.l 


2008; 


Modiaz et al.l 120081 


Levesaue et al. 12010a: 


Chornock et al.l 120101: IStarling et al.l 1201 ID. 


Chal- 



lenging this view is the recent discovery of SN2009bb, 
a broad-lined, engine-driven Type-Ic supernova found 
in a h igh-mcta llicity host env ironment (Sodcrberg et all 
120101 : [Levesa ue et al.l l2010dl : E^iiata et ah. 2011,) . In 
SN2010ay, we have found the opposite case - a broad- 
lined Type Ic supernova found in a low-metallicity host 
environment, but without any indication (via either 
radio or gamma-ray emission) of a central engine. 
The existence of SNe 2009bb and 2010ay emphasizes 
that progenitor metallicity is not the key factor that 
distinguishes GRB-SNe from broad-lined SNe Ic without 
associated relativistic ejecta. 

We compare the absorption velocity of SNe Ic-BL and 
engine-driven SNe (GRB-SNe and SN2009bb) to the 
metallicity of their host environments in Figure [TOl This 
comparison emphasizes the diversity of explosion and 
host galaxy properties observed in both engine-driven 
SNe and SNe Ic-BL. The engine-driven SNe with the 
largest velocity gradients (2003dh, a — —0.9 ± 0.1; 
2009bb, a = -0.86 ± 0.08) occur at metalhcitics dif- 
ferent by a factor of 5 (2003dh, log(0/H)-|-12 = 8.0; 



16 



Sanders ct al. 



# Ic-BL 

# GRB-SNe 



SN 2009bb 
SN 201 Oay 



SN 201 Obh 




-15 -16 -17 -18 -19 -20 -21 -22 -23 

Host Mb 



Fig. 9. — A plot of host galaxy metallicity versus absolute B magnitude for SNe Ic-BL (red) and engine-driven explosions (blue). The L-Z 
relation of nearby starforming galaxies is plotted as a solid line, v^ith the 15th and 85th percentile boundaries of the galaxy distribution 
(dash ed lines). Here wc have transformed the fc-corrccted gri m agnitudes from the MPA/JHU catalog to B-band i lBlanton Ik RoweisI 
120071) and converted the T04 metallicity values to the PP04 scale IIKewlev fc Ell ison 2008), for the purpose of comparing it to metallicity 
measurements for SN host galaxies in the literature. The dot-dashed horizontal line is the divider betwreen GRB-SNe and SNe Ic-BL 
host envir onments suggested b yl^Iodjaz ct al. (2008). The host galaxy prope rties of GRB/SNe other than 2 010ay arc fro m the following 
referen ces: IStarline et al.l2011l : [(Tano ct al..2011a f2010bh'). [T^odiaz et al.l2011l f Ic-BL). rCevesgue et al.l2010dl (2009bb), and lLevesQue et al.l 
I2010al (other GRB-SNe). Errorbars illustrate measurement uncertainty, when published, plus a 0.07 dex systematic uncertainty. 



2009bb, log(0/H)-fl2 = 8.7 on the PP04 scale). Fur- 
thermore, the velocity gradient of engine-driven SNe 
seems to be uncorrclated with the velocity at late times; 
2010bh {a = -0.22 ± 0.07) and 2006aj {a = -0.3 ± 0.2) 
have similar velocity gradients, but 2010bh had veloc- 
ities ~ 10, 000 km s~^ larger at 30 days after explo- 
sion. Among SNe Ic-BL from low metallicity environ- 
ments (log(0/H)+12 < 8.5), there is a large range in 
both the characteristic velocity and velocity gradient 
= [10 ± 1, 9 ± 2, 21 ± 2] and a = [-0.5 ± 0.2, -0.3 ± 
0.2, -0.4 ±0.3] for SNe [2003jd,2007bg,2010ay]). For the 
three SNe Ic-BL at higher metallicities, the character- 
istic velocities tend to be lower and the velocity gradi- 
ents tend to be stronger (v|P = [6 ± 2, 5 ± 2, 10.3 ± 0.7] 
and a = [-0.8 ± 0.4,-1.5 ± 0.7,-0.5 ± 0.1] for SNe 
[1997ef,2002ap,2007ru]). However, a larger sample is 
needed to exclude the possibility of SNe Ic-BL from 
super-solar metallicity environments that have high char- 
acteristic velocities or shallow velocity gradients. 

The fact that a GRB was not detected in association 
with SN2010ay, despite its similarities to the GRB-SNe, 
could indicate that the relativistic jet produced by this 
explosion was "s uffocated" before it em erged from the 
progenitor star (jMacFadven et aH l2001[ ). In this sce- 
nario, the duration of the accretion event onto the newly- 



formed central engine is short-lived and the jet post- 
breakout outflow is not ultra-relativistic. In the process 
of being suffocated, the jet transfers momentum to the 
ejeeta such that the spectrum is broad lined and the 
ejeeta velocity is very high, even at late times, as we 
observe (Q. However, the low host environment metal- 
licity we measure for SN 2010ay; which is similar to GRB- 
SNe, does not suggest high angular momentum loss that 
would help to weaken the jet. Another alternative is that 
binary interaction plays a key role in the commonality of 
relativistic outflows among SNe Ic-BL. 

Looking forward, additional SNe Ic-BL in sub-solar 
metallicity host environments will be found through 
current and future generations of untargeted transient 
searches. Untargeted searches arc required to find su- 
pernovae from low-metallicity host environments, be- 
cause targeted searches only probe the highest metallic- 
ity galaxies due to the lum inosity-metallicity relationship 
(jModjaz et al.ll2011l : iLelo udas et al. 2011,). 

8. CONCLUSIONS 

The optical photometric and spectroscopic, radio, and 
gamma-ray observations of SN 2010ay presented here 
provide an example of a Type Ic-BL SN with explo- 
sion and host properties similar to the known GRB-SN 
SN2010bh. This object demonstrates that SNe in low- 
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T 




log(0/H+12) (PP04 03N2) 



Fig. 10. — SN absorption velocity, as traced by the Si II A6355A feature, versus host galaxy oxygen abundance for SN2010ay and other 
Ic-BL (red) and engine-driven explosions (blue) from the literature (as in Figure [Sjl. The range of velocities hatched for each object comes 
from the velocity at 10 days and at 30 days after explosi on, according to the po wer law fits p resented in Figure [S] The oxyge n abundance 
measu rements using the PP04 03N2 d i agnos tic are from ILevesgue et al.l Il2010al ) (GRB-SNe). ISahu et all 120091) (SN2007ru'). lYoung etaU 
1)20101 ) (SN 2007bg), and lModiaz et al] l|20111 ) (other SNe Ic-BL). The range of oxygen abundanc e hatched reflects the er ror bars quoted in 
the literature (when stated) plus the ~ 0.07 dex systematic error of the PP04 03N2 diagnostic IjKewlev Ik Ellisonll2008l ). 



metallicity environments with high-velocity ejecta are 
not necessarily accompanied by the traditional signa- 
ture of radio emission associated with long-lived relativis- 
tic jets. The existence of SN 2010ay and SN 2009bb (a 
central-engine driven event from a high-metallicity host 
environment) indicates that progenitor metallicity may 
not be the key factor that distinguishes GRB-SNe from 
normal broad-lined Type Ic supernovae. 
We conclude that: 

1. Pre-discovery imaging of the SN2010ay from the 
Pan-STARRSl Stt survey allows us to tightly 
constrain the early-time light-curve of SN 2010ay 
(see Figure [5]) and explosion date (2010 February 
21.3±1.3), allowing us to search for gamma-ray 
emission that may have been associated with the 
explosion. By fitting the template SN Ib/c light- 



curve of iDrout et all (poTTI ). we derive an R-band 
peak absolute magnitude of — 20.2±0.2 mag - mak- 
ing SN2010ay among the most luminous SNe Ib/c 
ever observed. This peak magnitude suggests that 
a large mass of nickel, il/Ni ~ 0.9^0 ;[ 
synthesized. We estimate a ratio of A/n 
is '--^ 2x larger than in known GRB-SNe. 



0.9Ip;i Mq, has been 
to il/oj that 



2. Spectroscopy (see Figure ^ at the explosion site 
in the host galaxy of SN2010ay indicates that the 
host environment of the progenitor star had a sig- 
nificantly sub-solar metallicity (Z 0.3 Zq), sim- 
ilar to the host environments of known GRB-SNe 
progenitors. 

3. The Type Ic-BL SN 2010ay strongly resembles the 
GRB-SN2010bh, particularly in light of its un- 
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usually high absorption velocities at late times 
{vsi «19.2xl0^ km s~^ at 14 days after peak) and 
low metallicity host environment. The comparison 
between these two SNe is summarized in Table [S] 

4. Non-detections in late-time EVLA radio observa- 
tions of the SN rule out the association of a GRB 
of the nature of the spectroscopically-confirmed 
GRB-SNc. except for the radio afterglow associ- 
ated with XRF 060218. Our radio observations im- 
ply limits on the velocity, energy, and density of any 
associated relativistic jet and the mass loss rate of 
the progenitor (see Figures HI El [71 and [8]) . Ad- 
ditionally, no coincident gamma-ray emission was 
detected by satellites: the non-detection by the in- 
terplanetary network indicates < 6 x erg, 
while the non-detection by the Swift BAT indi- 
cates that the peak energy of the burst was < 
1 X 10"*^ erg s^^ if the burst occurred during the 
~ 20% of the explosion window when it was in the 
field of view of the instrument. This rules out as- 
sociated prompt emission similar to that of GRBs 
031203, 030329, or 100316D, but not GRBs 980425 
or 060218. 

The pre-discovery imaging of SN 2010ay demonstrates 
the capability of the untargeted PSl survey for identi- 
fying and monitoring exotic transients, not only in its 
high-cadence Medium-Deep Fields, but also in the all-sky 
3tt survey. Additional detections and multi-wavelength 
follow-up observations of SNe Ic-BL will help to illumi- 
nate the role that ejecta velocity and progenitor metal- 
licity play in the GRB-SNc connection. 
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